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Abstract 

A method is proposed for estimating the potential function of a 
non-parametric estimator for stationary and isotropic pairwise inter¬ 
action point process. The relation between a pair potential and the 
corresponding Papangelou conditional intensity is considered. Consis¬ 
tency and strong consistency of non-parametric estimate are proved in 
case of finite-range interaction potential. 
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1 Introduction 


Gibbs point processes are a natural class of models for point patterns exhibit¬ 
ing interactions between the points. By far the most widely applied form in 
practical analysis is that of pairwise interaction, where the scale and strength 
of interaction between two points are determined by a so-called pair potential 
function. For a stationary and isotropic process the pair potential is a func¬ 
tion of the distance between the two points. Fields of applications for point 
processes are image processing, analysis of the structure of tissues in medi¬ 
cal sciences, forestry (Matern [?]), ecology (Diggle [?]), spatial epidemiology 
(Lawson [?])and astrophysics (Neyman and Scott [?]). 

Pairwise interaction point process densities are intractable as the nor¬ 
malizing constant is unknown and/or extremely complicated to approximate. 
However, we can resort to estimates of parameters using the conditional in¬ 
tensity. In this paper, we suggest a new non-parametric estimate of the 
pair potential function for stationary and isotropic pairwise interaction point 
process specihed by a Papangelou conditional intensity on increasing regions 
single realization is observed. In this cas a point process is dehned as a ran¬ 
dom locally-hnite counting measure on the d-dimensional Euclidean space 

. Consistency and strong consistency of the resulting estimator are es¬ 
tablished. 

To our knowledge only one attempt to solve the problem of non-parametric 
estimation of the pair correlation function and its approximate relation to 
the pair potential through the Percus Yevick equation (Diggle et ah [?]). 
The approximation is a result of a cluster expansion method, and it is ac¬ 
curate only for sparse data. Many attempts have been tried to estimate 
the potential function from point pattern data in a parametric framework ; 
maximization of likelihood approximations (Ogata and Tanemura [?], Ogata 
and Tanemura [?], Penttinen [?]), pseudolikelihood maximization (Besag et 
ah [?], Jensen and Mpller [?]) and also some ad hoc methods (Strauss [?], 
Ripley [?], Hanisch and Stoyan [?], Diggle and Gratton [?], Fiksel [?], Takacs 
[?], Billiot and Goulard [?]). 

Our paper is organized as follows. Section |2] introduces basic notation 
and definitions. In Section [3l we briefly present some models satisfying the 
assumptions needed to prove our asymptotic results. In Section 01 we present 
our main results. Consistency of non-parametric estimator is proved in Sec¬ 
tion 01 it is based on the knowledge of Papangelou conditional intensity and 
the iterated Georgii-Nguyen-Zessin formula. Using Orlicz spaces we can ob¬ 
tain a strong consistency of non-parametric estimator in Section [6l 


2 


2 Basic notation and definitions 


Throughout the paper we adopt the following notation. We denote the space 
of locally hnite point conhgurations in by Nif. The volume of a bounded 
Borel set W of R'^ is denoted by |hT| and o denotes the origin. For all hnite 
subset T of we denote |r| the number of elements in T. || ■ || denotes 
Euclidean distance on R'^. ad = r(d/ 2 ) measure of the unit sphere in 

R'^. Let be the unit sphere in R^. 

Papangelou conditional intensity (Mpller and Waagepetersen [?]) of pair¬ 
wise interaction point process has the form 

A(u,x) = 7 o(M)exp ( - 7 o({u,u}) 

^ vG'x\u 

If 7 o(u) = /9 is a constant and 7 o({m,u}) = 7 (||u — u||) is invariant under 
translations and rotations, then a pairwise interaction point process is said 
to be stationary and isotropic or homogeneous. The Papangelou conditional 
intensity can be interpreted as follows: for any u G R'^ and x G Nif, X{u, x)dM 
corresponds to the conditional probability of observing a point in a ball of 
volume du around u given the rest of the point process is x. Fortunately 
does not contain a normalising factor. 

For convenience, throughout in this paper, we consider stationary and 
isotropic pairwise interaction point process. Then its Papangelou conditional 
intensity at a location u is given by 

AKx)=/3*exp('- 7 (||t;-«||)Y VugR",xgW/ (2.1) 

where P* is the true value of the Poisson intensity parameter , 7 is called the 
pair potential, a name that originates in physics: it measures the potential 
energy caused by the interaction among pairs of points {u, v) as a function 
of their distance ||u —u||. Usually a hnite range of interaction, R, is assumed 
such that 

7 (||u —m||) = 0 whenever I|u — u|I > i?. ( 2 . 2 ) 

We assume that 7 (||u — u||) > 0 for ||u — u|| < R, so that typical realizations 
will be more or less regular compared to a completely random arrangement. 
The pairwise interaction between points may also be described in terms of 
the pair potential function 7 into the interaction function <P = exp(— 7 ). For 
> 1, A(u,x) is increasing in x. For ^ < 1, A(m,x) is decreasing in x (the 
repulsive case). It can be computed for the case ^ = 1 which corresponds to 
the homogeneous Poisson point process with with intensity /9*. 
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3 Examples of Papangelou conditional inten¬ 
sity 

Examples of conditional intensities are presented in Baddeley et al [?], Mpller 
and Waagepetersen ([?], [?]). The following presents some examples which 
have been applied in various contexts and satisfying the assumptions needed 
to prove our asymptotic results. 

1. A special case of pairwise interaction is the Strauss process. It has 
Papangelou conditional intensity 

A(n,x) = 

where /9>0,0<^<1 and n[o,_R](M,x) = ^Il(||n —n|| < R) is the 
number of pairs in x with distance not greater than R. 

2. Piecewise Strauss point process. 

P 

A(n,x) = /3j_p. ^ ^ 

.7 = 1 

where P > 0, 0 < < 1, x) = ^ l(|ln - u\\ e [Rj-i, Rj]) 

v£x. 

and Rq = 0 < Ri < ... < Rp = R < oo. 

3. Triplets point process. 

A(m,x) = /5^*[0,fl](^U77)-S[0,H](xV) 

where /9>0, 0<^<1 and S[o,/j](x) is the number of unordered 
triplets that are closer than R. 

4. Lennard-Jones model 

A(m,x) = /3 Yl ^{\\v-u\\) 

vGx.\u 

with log^(r) = for r = ||n —m||, where 6 > 0 

and P > 0 ares parameters. 
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4 Main results 


Suppose that a single realization x of a point process X is observed in a 
bounded window Wn where {Wn)n>i is a sequence of cubes growing up to 
Throughout in this paper, h is a non-negative measurable function dehned 
for all u G x G by 


note that 

and 


h{u, x) = 1 ("inf I |n — m| I > i?) =1 {d{u, x) > i?), 

\v£x / 


F{o, rv) = E[h(o, yi)h{rv, x) 


J{r)= / F(o, rn)dn. 


To estimate the function J(r)^(r), we introduce edge-corrected kernel- 
type estimator Rn{r) dehned by 


Rnir) = 


bn\WnQ2R\o-d ^^\\V - U 
u,v£X. 

Ilf— 


X] v})h{v, X\{n, v})Ki f 


— r 


(4.3) 


© will denote Minkowski substruction, with the convention that 

WnQ 2 R = hhn © B{u, 2R) = {u G Wn ■ ||n — n|| < 2R for all v G Wn} 
denotes the 2i?-interior of the cubes Wn, with Lebesgue measure \ WnQ 2 R\ > 0. 

signihes summation over distinct pairs, i^i : R —)■ M is an univariate 
kernel function associated with a sequence (&n)n>i of bandwidths satisfying 
the following: 


Condition i^(l,Q;) ; The sequence of bandwidths > 0 for n > 1, is 
chosen such that 


lim = 0 and lim hn\Wne 2 R\ = cxd. 

n^oo n^oo 

The kernel function Ki : R —)■ R is non-negative and bounded with bounded 
support, such that: 


Ki{u)du = 1, 


u^Ki{u)du = 0,j = 0,1,..., 0 — 1, for a > 2. 
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To estimate the function we introduce empiric estimator Jn{r) defined 

by 


^ - 1 X \ {u})h%u, X \ {«}), (4.4) 

WnQ2R\ 

where h*{u,x.) = f^d-i h{rv — M,x)dn. Using the spatial ergodic theorem of 
Nguyen and Zessin [?], estimator fl4.4p turn out to be unbiased and strongly 
consistent. The natural estimator of Poisson intensity (3* is 

_ E^6X 

r ut * V*/ 

This estimator turn out to be unbiased and strongly consistent and results 
on asymptotic normality were obtained by Morsli et ah [?]. 

Plugging in the above estimator and then the interaction 

function ^(r) = exp(— 7 (r)) for r G (0, i?] can be estimated using edge- 
corrected non-parametric estimate by 


Rn{r) 


(3nJn{X ) 


(4.6) 


The strong consistency of the estimators (14.41) and f|4.5p implies the following; 

Proposition 1. Let 7 he pairwise interaction potential defined in (12.ip sat¬ 
isfying condition ( 1 ^ . Let Ki kernel function satisfying Condition Kfi^a) 
and the function J{r)exp{—^{r)) has bounded and continuous partial deriva¬ 
tives of order a for all a > 1 in {r — S,r -\- 5) for some 5 > 0. Then as 
n ^ 00 


—texp(— 7 (r)) in probability P (respi¥-a.s.) iff 
2 

Rn{f) —J(r) exp(— 7 (r)) in probability P (respJS’- a.s.). 

The convergence in probability (consistency) for a wide class of point 
process will be discussed in Section [5l Conditions ensuring uniform P-a.s. 
convergence of kernel-type estimator of i?„(r) and the strong consistency 
<Ln{r) will be discussed in Section O 
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5 Consistoncy 


5.1 Asymptotic behaviour mean squared error of the 
kernel-type estimator 

In this section we will derive bounds for the mean squared error of the ker¬ 
nel estimator kernel-type estimator of Rn{r). We consider the mean square 
error of Rnir)^ MSE (i?„(r)) = Var (i?„(r)) -|- {^Biais{Rn{r)'f. So con¬ 
vergence in MSE implies that as n ^ oo [Biais{Rn{r)Y = (Ei?„(r) — 
/3*‘^J{r) exp(— 7 (r)))^ —)■ 0 and Var (i?„(r)) = E (-R„(r) — Ei?„(r))^ —)■ 0. 
Hence, -R„(r) is consistent in the quadratic mean and hence consistent. For 
doing this, we first determine the asymptotic behaviour of Ei?,^(r) and 
Var Rn{r). 

Theorem 1. Let 7 be pairwise interaction potential defined in fl2.ip satisfying 
condition fl 2 . 2 p . Let Ki kernel function satisfying Condition K{1,1). For all 
r G (0, R], we have 

lim ERn{r) = J(r) exp(— 7 (r)). 

n—>-cxD 

If Condition K{1, a) is satisfied and the function exp(— 7 (r)) J(r) has bounded 
and continuous partial derivatives of order a in {r — S,r + 6) for some 5 > 0 
and for all a >1. Then 

ERn{r) = l3*‘^J{r) exp(— 7 (r)) -|- as n ^ 00 . 

Theorem 2. Let 7 be pairwise interaction potential defined in fl2.ip satisfying 
condition ( 123 ). Let Ki kernel function satisfying Condition K{l,a) for all 
a > 1 such that Kf(p)dp < 00 . For all r G (0, R], we have, 

lim hn\Wne 2 R\^sx{Rn{r)) = J(r)exp(- 7 (r)) [ K^{p)dp. 

n^oo CTdr^-^ Jr 


5.2 Proof of Theorem [T] 

Proof. We define 

L{ui,..., Us, X) = h{ui, X)...h{us, X), F{ui,..., u^) = E[h{ui, X)---h{us, X)] 
and J(||mi||, •••, ||m^||) = 1(||mi|| < R, •••, ||w|| < R)- 
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The calculation of expectation and variance of -Rn(r) is based on the 
iterated Georgii-Nguyen-Zessin (GNZ) formula, see Papangelou [?]; 


h{ui, {ui, 

ni,...,UsGX 


E h{ui, tts, X)A(ni,tts, X)dtti...dn„ 

(5.7) 


for non-negative functions h : x Nif —> M, where A(mi, x) is 

Papangelou conditional intensity and is dehned (not uniquely) by 

A(mi, ...,Ms,x) = A(Mi,x)A(n2,xU {mi})...A(m5,xU {ui, 

Applying the preceding formula fl5.7p for s = 2, we derive 

^Rn{r) = I / -^ 

0n\Wne2R\o-d 


E 




\v — u\ 


^J(||n-n||)L(n,n,X)iri 


\v — u\ \ — r 


A(m, V, X)dMdn. 


For an interaction radius R, the Papangelou conditional intensity satisfies 

A(m,x) = A(m, 0) for all x with d{u,x.) > R 

since points further away from u than R do not contribute to the Papangelou 
conditional intensity at u. Using the hnite range property fl2.2|) . we get 




*2 


RRnir) = - 

&n|hhne2ij|c’‘d 


E 




I 

3*2 


\v — u\ 


r^Ji\\v-u\\)L{u,v,X)K, 




bn 


A®E|L(o,s,X)1A-, 


\v — u\ \ — r 


^(||s||)ds. 


^{\\v — nlDdudn 


Recall a property of the integration theory (see Briane and Pages [?] or 
Rudin [?]). Let be the unit sphere in i.e. E>^~^ = {u G : ||m|| = 1}, 
then for any Borel function / : — 




f{u)du = 


f{rz)r'^ ^addrdz. 


'0 Jsd-^ 

By combining the above result, we get so: 

/ 5*2 p . 


ER„(r) = 


J-r/bn 


J{hnQ + r)F(o, {hnQ + r)v)Ki{Q)${hnQ + r)addQdv. 














With bounded support on the kernel function and by dominated convergence 
theorem, we get as n —)■ oo, Ei?„(r) — > J(r) exp(— 7 (r)). Now, we are 
going to prove the second part of the Theorem [T] We have a product of two 
functions F{o, {hnQ + r)v)<P{hnQ + t) and we approximate each one of them 
with a Taylor formula up to a certain a. We use Taylor’s formula to obtain 
for n ^ oo, 


a—1 


<P{bng + r) = ${r) + ^ 


(bng)^ 

^r)dr + ^^-—{r + hr,Qe) 


k=l 


k\ dr 


al dr'^ 


and 


Q —1 


F{o, {bnQ + r)v) = F{o, rv) + ^ 


{bng)^ dF{o,rv) b^d°‘F 


k=l 


k\ 


dr 


dr + 


a\ dr'^ 


(o, (r + bng9)v). 


So we denote this product by T^irv, r), then we have as n —)■ oo 


a—1 

F{o, [b^g + r)v)<P{bng + r) = F{o,rv)<P{r) + ^ T„(ru, r)(6np)^ + 0(&“). 

k=l 


It follows that, 

E^„(r) = /3*V(r)^(r) 


a—1 


+ / y'^n(ru,r)6^du 

+ C>(6“) as n —)■ oo. 


g'"Ki{g)dg 


Together with Condition i^(l,Q!) imply the second assertion of Theorem 

m ■ 


5.3 Proof of Theorem [2] 

Proof. The proof of Theorem [2] makes use of the following corollary. 

Corollary 1. Consider any Gibbs point process X in with Papangelou 
conditional intensity A. For any non-negative, measurable and symmetric 
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function f : Mf' x x Nif — > M, we have 
Var^ f{u,v,X.\{u,v}) 

^u,vex 

= 2E f f‘^{u,v,X.)X{u,v,X.)dudv 

jR2d 

+ 4E / f{u,v,X.)f{v,w,X.)X{u,v,w,'K.)dudvdw 

Jm.3d 

+ E / f{u,v,X.)f{w,y,X.)X{u,v,w,y,X.)dudvdwdy 
- [ E[/(M,i;,X)A(M,i;,X)] E[/(w;,|/,X)A(w;,l/,X)]dMdwdw;d|/. 

AR4d 

Proof. Consider the decomposition (see Jolivet [?] and Heinrich [?]) 

7^ \ 2 ^ 

Y j =2 ^ /^(M,'i;,X\{M,r;}) 

+ 4 5]] /(M,n,X\{M,n,M;})/(n,M;,X\{M,n,M;}) 

w,i;,tL!GX 

7^ 

u,v,w,y^'X. 

(5.8) 

Applying the preceding (GNZ) formula fl5.7p combining with fl5.8p . we obtain 

7^ 

Var Y fiu,v,X\{u,v}) 

ti,wGX 

7^ 7^ 

= E( “ (E Y 

hijh’GX u,vG'K 

= 2E f f‘^{u,v,X)X{u,v,X)dudv 
JR2'^ 

+ 4E / /(m, n,X)/(n,tc,X)A( m, n,w,X)dMdndic 

jR3d 

+ E / /(m, n, X)/(r(;, ?/, X)A(m, n, ta, ?/, X)dMdndr(;d|/ 


E[/(m, n, X)A(m, n, X)] E[/(u;, ?/, X)A(m;, ?/, X)]dMdndu;d|/. 
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We obtain the desired result. 


Applying Corollary [T] to this function 


f{u,v,X) = 




\v — u\ 


j^J{\\v-u\\)L{u,v,X)K^ 


\v — u\ 


bime2R?aj 


it is easily seen that Vari?„(r) = Ai + A 2 + As — A^, where 

E 

A 2 = 


lWne2H(“) 


\v — u 


''^^^''^7^J{\\v — u\\)L{u,v,X)K‘^[ b--) A(n, n, X)dMdn, 


||2(d-l)' 


\v — u\ 




WW7 -/ II I,; /W (||n-n||,||n-n;||)LKn,n;,X) 

bl\Wne 2 R\^(Tj U, \\v - uW^-^v - w\\^-^ 


Ki 


— 


\v — u\ 


Ki 


W — w\ 


■E 


bl\WnQ 2 R\'^(Tj jRid lb - u\\<^~mw - y\\ 


X{u, V, w, X)dndndta, 


^wyi,e2fl(^)^n^n,e2fl(^) T/|| II II \\\Tf 

J(|b- m||, ||w- 2 /||)L(M,n,n;,|/,X) 


d-l 


X ATi 

and 
A 4 = 


\v — u\ 


Ks 


\w 


\{u, V, w, y, X)dMdndtcd|/, 


bime2RK 


E 


lvyne2H(“) 


Q 2 R y 

\v — u\ 


^J{\\v-u\\)L{u,v,X)Ks 


\v — u\ \ — r 


X{u, V, X)dMdn 


V. 




The asymptotic behaviour of the leading term Ai is obtained by applying 
the second order Papangelou conditional intensity given by: 

A(m, n, x) = A(m, x)A(n, X U {«}) for any n, n G and x G W/- 

Using the finite range property fl2.2p for each function A(m, x) and A(n, xU 
{«}), this implies that 

A(m,0) = /9* and A(n, 0 U {«}) = /3*^(| jn — m| |) for all M,nGM‘^. 
And by stationarity of X, it results 
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Ai — 


2/3' 


'*2 


bl\W^e2R?aj 




»»»'“) j(||„ _ .„||) EIH{o,v- u,X)]K'f 


l^ 2 d \ \v - 
2/3*2 


\v — u\ \ — r 


^(llu — u\\)dudv 


bn\Wne2R\cTa J_r/b„ {buQ + r) 


+ + r)d^dcr(M;). 


Dominated convergence theorem and assumption of Ki imply for all r G 

{0,R] 

lim bn\Wne 2 R\Ai = J(r)^(r) [ K‘^{p)dp. 

n^oo CTrfr‘^-1 Jjj 

We will now show that all other integrals to Var i?„(r) converge to zero. 
For the asymptotic behaviour of the second term A 2 , we remember the third 
order Papangelou conditional intensity by 


X{u,v,w,x.) = A(m, x)A(u, X U {m})A(w, x U {u, u}) 


for any u,v,w G and x G iV;/- Since X is a point process to interact 
in pairs, the interaction terms due to triplets or higher order are equal to 
one, i.e. the potential 7 (y) = 0 when u(y) > 3, for 0 7 ^ y C x. Using 
the finite range property fl 2 . 2 l) for each function A(m,x), A(u,xU {m}) and 
X{w,:k U {m,u}) and after a elementary calculation, we have 


X{u,v,w, 0) 


/3*^(p/\\v — u\\)${\\w — v\\) if d{u,w) < R 
/3*^^(||u — m||) otherwise. 


Which ensures that X{u,v,w,(/)) is a function that depends only variables 
||u — M||,||tc — u||, denoted by ^i(||u — m||, ||tc — u||). 

According to the stationarity of X, it follows that 


Ao — 




IFTFiT-/ II ill I / 

bl\Wne2R\^(Tj Ud \\v - uW^-^v - w\\^-^ m; V , , , 


X ^i{\\v — u\\,\\w — v\\)KI 
4 


\v — u\ 


Ki 


\v — w\ 


■) 


dududtc 


^CXD poo 


\WnQ2R n {Wne2R “ [b^Q + r)z)\ 


\WnQ2RW‘d J-r/b„ J-r/b„ |fKie 2 _R| 

X F{o, {bnQ + r)z, {bng' + r)z')Fi{bng + r, + r)Ki{g)Ki{g')dgdg'dad{z)dad{z'). 
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The asymptotic behaviour of the leading term A 2 is obtained by applying 
the dominated convergence theorem. When multiplied by 6 „lTne 2 R|, we get 
\\m.n^oohn\Wne2R\A2 = 0. 

Next we introduce the finite range property f|2.2|) and reasoning analogous 
with the foregoing on X{u,v,w,y,(/)). which ensures that X{u,v,w,y,(/)) is a 
function that depends only variables | |n — nl|, 11|/ — tc||, ||w — nl |, | |w — n| |), 
denoted by ^ 2 {.\\v — nil, \ \y — tc||, ||tc — m||, ||tc — n| |). We find that 


A f 17 /" ('^) j .^11 II II ii\r/' 

^3 = TyrfTT--r[^d( n - t 6 , - y\\)L[u,v,w,y,X) 

bi\Wne2Rr(rj J^-i\\v-u\\^ - y\\‘^ i 


X Ki 


V — u — r 


Ki 


\w — y\ \ — r 


\Wne2RWj 


[ r f I i 

jRd J-rlbn J-r/bri 


X{u, V, w, y, X)dndndr(;d?/ 
\WnQ2R n {WnQ2R “ 'n>)| 


\WnQ2R\ 


Ki{q)K^{q') 


X <^2{bnQ + r, bnd + r, ||n;||, ||(6„£i + r)z - w\\)dwdQdQ'dadiz)dadiz'). 


Where 

^* 2 {bn 0 +r,br,g' + r,\\w\\,\\{bng+r)z-w\\) = J(6„p+r, 6„p' + r)F(o, (6n^+ 
r)z, {bng' + r)z')^2{bng + r, bng' + r,\\w\\,\\{bng + r)z - w\\). 

By dominated convergence theorem, we get lim„_j.oo&n|hfne 2 i?l ^3 = 0- 
For asymptotic behaviour of the leading term ^ 4 , it then suffices to repeat 
the arguments developed previously to conclude the following result. 


Aa = 




*4 




[ 


X)1 E|i(», y. X)1 


X ^{\\v - n||)^(||y - n;||)ifi 


\v — u\ \ — r 


Ki 


\w — y\ \ — r 


, I dndndrcd?/ 

bn J 


[ 


bl\WnQ 2 R\(Tj 

|lhne 2 j; n {WnQ2R — w)| 
\Wne2R\\\v - u\\<^-^\\w - 


^ J(||n||, I In; - y||) E[//(o, n, X)] E[L(n;, y, X)] 


X ^>(||n||)^(||y - u;||)Xi 


m — r 


Ki 


|n; — y|| — r 


dvdwdy 


\Wne2R\(Tj 

r-- roo r r \Wne2R n iWne2R - w)\ 


r r r r r 

J —rlhn J —rfb-n ^ ^ 


\Wne2R\ 


J {bnQ + r, bnd + r)F{o, {bnQ + r)z) 


X Fo, {bng' + r)z')<P{bnQ + r)F{bng' + r)Ki{g)Ki{g')dwdgdg'dad{z)dad{z'). 
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Then by dominated convergence theorem, we get hm„^oo &n|hh’„e 2 R|A 4 

0 . 


6 Strong consistency 

6.1 Rates uniform strong convergence of the kernel-type 
estimator 

Before realizing the strong consistency ^n(^) we introduce some necessary 
dehnitions and notation. A Young function is a real convex nondecreasing 
function dehned on which satishes hmj_j,oo = +oo and '^(0) = 0. We 
define the Orlicz space as the space of real random variables Z dehned 
on the probability space (W/, P) such that E[ip{\Z\/c)] < +cx) for some 
c > 0. The Orlicz space equipped with the so-called Luxemburg norm 
||.||^ dehned for any real random variable Z by 

\\Z\\^ = mi{c>0;E[ij{\Z\/c)]<l} 

is a Banach space. For more about Young functions and Orlicz spaces one 
can refer to Krasnosel’skii and Rutickii [?]. Let 6 > 0. We denote by ipQ the 
Young function dehned for any x G M"*" by 

'ipeix) = exp((a; - exp(^^) where ie = ((1 - <6 <1}. 

On the lattice we dehne the lexicographic order as follows: if f = (ii,..., 
and j = {ji, are distinct elements of the notation i <iex j means 

that either R < ji or for some p in {2, 3,..., d}, ip < jp and iq = jq for 
1 <q < p. Let the sets {V^ ] i E , k eN*} be dehned as follows: 

Y/ = {j e Z^ J <iex i}, 


and for k >2 

n {j e Z'^; \i — j\ > k} where \i — j\ = max \ii — ji\. 

l<l<d 

For any subset F of Z*^ dehne Er = a(ei; f G F) and set 

E\k\{£i) = E{ei\Ey\k\), k E . 

Denote 9{q) = 2g/(2 — q) for 0 < g < 2 and by convention l/9{2) = 0. 
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Next we list a set of conditions which are needed to obtain (rates of) 
uniform strong consistency over some compact set [ri,r 2 ] in ( 0 , i?] of the 
estimator -R„(r) to the function J(r)^(r). The following assumption is 
imposed; 

Condition Cp : The kernel function is a Lipschitz condition, i.e. there 
exists a constant 77 > 0 such that 

\Ki{p) - Ki{p')\ <r]\p-p'\ for all p, p'e [ri,r 2 ]. 

Strong uniform consistency for the resulting estimator are obtained via 
assumptions of belonging to Orlicz spaces induced by exponential Young 
functions for stationary real random fields which allows us to derive the 
Kahane-Khintchine inequalities by El Machkouri [?]. Our results also carry 
through the most important particular case of Orlicz spaces random fields, 
we use the inequality follows from a Marcinkiewicz-Zygmund type inequality 
by Dedecker [?]. 

Now, we split up the sampling window WnQ 2 R into cubes such as WnQ 2 R = 
UieFnAj, where Aj are centered at i and assume that T^ = {—n, ..., 0 , 
increases towards Z*^. We split up Rnir) as follows: 

Rn{r) = ^— 'V] Rk{r) 

bn\Wne2RWd ^ 


U \ Oyi } 

«,t>ex 11 11 \ n / 

\\v—u\\<R. 

Note for all k G T^, Rk = Rk{r) - RRk{r) and Sn = Z]fcGr„ Rk{r). 

Theorems. Under Conditions K{1, a) andCp. Further, assume that J{r) exp{—'y{r)) 
has bounded and continuous partial derivatives of order a in [ri — 6, r 2 + 5] 
for some <5 > 0 . 


1) If there exists 0 < q < 2 such that Rq G and 


E 

kevf 


RkR\k\ (-^ 0 ) 


< CX). 


'le(q) 


(6.9) 


Then 

sup |.R„(r)-/3*V(r)exp(-7(r))| = C>a.s. )+C»(&n) «« nex. 

ri<r<r2 \ ' / 
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2 ) If Rq G L°° and 


< oo. 


( 6 . 10 ) 


Then 


\\RkE\k\{Ro) 


sup |.R„(r)-/?*V(r)exp(- 7 (r)) 

ri<r<r2 


/ {fognyR 

\ bnn‘^1'^ 


+C?(&n) 


as 


n ^ oo. 


3) If there exists p > 2 such that Rq G and 

\\RkE\k\{Ro)\\E < oo. (6.11) 

Assume that bn = n“'^^(logn)'^i for some qi, q 2 > 0- Let a,b > 0 be fixed 
and if a{p + 1) — df/2 — q 2 > I and b{p + 1) + gi > 1. Then 

sup [.Rn(r)-/3*V(r) exp(- 7 (r))| = Oa.s.i ^ ^+0{b°) as n ^ oo. 

Remark 1. From the Markov property of X entails that for i ^ 0 are not 
neighborhoods, then Ri et Ro are conditionally independent, i.e'Ej[RQ\{X\r,i ^ 

0] = 0. Since a{Ri,i G Vq) is contained in a(XA.,i ^ 0) for k> I, for some 
integer I, it follows immediately that conditions fl6.9p . fl 6 . 10 D . fl 6 . 1 ip are sat¬ 
isfied. 


6.2 Proof of Theorem [3] 

Proof. To establish rates of the uniform P- a.s. convergence for the estimator 
Rn{r), we apply a triangle inequality decomposition allows for 




_ 

sup 

Rn{r) - ERn{r)\ < sup 

Rn{r) - Rn{p) 

Si-i<r<Si 

Si-i<r,p<Si 



+ 

+ 


sup 

Si-i<r,p<Si 


sup 

Si-l<p<Si 


ERn{r) -ERn{p) 

Rn{p) -ERn{p) 


The compact set [ri,r 2 ] is covered by the intervals Ci = [si_i — s*], where 
Si = ri i{r 2 — ri)/N,i = 1,...,N. Choosing N as the largest integer 
satisfying N < c/ln and Under the condition Cp, we deduce that 

there exists a constant 7 > 0 such that for any n sufficiently large 
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sup 

Si-i<r,p<Si 


Rn{r) - Rn{p) 


< -^'n\r- p\Rr, 

_ ' nA ^T). 


where 


7^ 


Rn = 




_^y 

<7j|»'„e2Hl IK’-« 

||ij—n||<_R 


T^h{u, X \ {m, v})h{v, X \ {m, u}). 


Follows from the last inequalities and the Nguyen and Zessin ergodic theorem 


sup 

Si-i<r,p<Si 


Rn{r) - Rn{p) 


= Op^si'^^'n) as n —)■ cxD. 


As well 


sup 

Si-i<r,p<Si 


ERn{r) -ERn{p) 


= Op,s.('>^n) as n —)■ cxo. 


Lemma 1. Assume that either O holds for some 0 < q < 2 such that 
Rq G and = (logn)^/'^/6„(y^)'^ or flh.lOp holds such that Rq G Loo 

and Un = {log/bn{^/nY. Then 


sup 

Si-l<P<Si 


Rn{p) - ERn{p) 


= Op,s.{rn) as n —)■ cxD. 


Proof. For e > 0, using Markov’s inequality, we get 


Rnir) -ERn{r)\ > er„ ) = P (|S„| > ernbnW 


< exp 


ernbnU^ 


ll^nllp, 


+ iq 


'e(i) 


Eexp 


\Sn 


WSnWp, 


+ 


SM 


Therefore, we assume that there exists a real 0 < g < 2, such that 
Rq G and using Kahane-Khintchine inequalities (cf. El Machkouri [?], 

Theorem 1), we have 


P ( \Rn{r) - Ei2„(r)| > er„ ) < P {^\Sn\ > ernbnn"^ 

< (1 + 6xp 


ernbnU^ 


b^AR)y/^ 


+ Co 
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denote 


k,iR)= poii;,„+E 




RkE\k\{Ro) 


i’e(q) 


We derive that if condition 06.91) holds, then there exist constant C > 0 and 
so if Tn = {logny/'^/bniyRlY 


sup P(|i?„(r) — ERn{r)\ > er„) < (1 + exp 

ri<r<r2 


e'^ logn 
Cl 


Now, we will accomplish the second step the proof of Proposition [3l Using 
Kahane-Khintchine inequalities (cf. El Machkouri [?], Theorem 1) with q = 
2, such that Rq G Lqo, we have 


P \Rn{r) - ERn{r)\ > ernj <2 exp 


£ Tnbnn'^ 


V(E«r. 


denote 

bi,2{R) = ||-Rol|^+ ^ \\RkE\k\iRo)\\^- 

fcGVol 

We derive that if condition Oh.lOp holds and so if = {lognYR/bn{y/nY, 
there exists C > 0 such that 


sup P(|.R„(r) - E.R„(r)| > er„) < 2 exp 

ri<r<r2 


logn 


choosing e sufficiently large, therefore, it follows with Borel-Cantelli’s lemma 


P(lim sup sup 

n—^■oo Si— 


Rnip) 


ERnip) 


> eVn) = 0. 


Now, we will accomplish the last step the proof of Theorem |3l 

Lemma 2. Assume flh.lip holds for some p > 2 such that Rq G and 
bn = n“® (logn)'^i for some constants qi,q 2 > 0. Let Let a,b > 0 be fixed and 
denote = n“(logn)^/6„(i/n)'^. If 

a{p + 1) - d/2 - q 2 > I et b{p + 1) + qi > 1, 


then 


sup 

Si-l<p<Si 


Rnip) 


ERnip) 


Cp.sivn) as n —)■ oo. 
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Proof. Let p > 2 be fixed, such that Rq G and for any £ > 0, 

P(|^„(r) - E^„(r)| > evn) = P (|S'„| > ernhnu'^) 

e-^R\Sn\^ 


~ rWnUP^ 
T-P 


— .r-P 


e 


rnb^nP'^ 


‘^V^Ci{R) 


p/2 


ier„ 


The last inequality follows from a Marcinkiewicz-Zygmund type inequality 
by Dedecker [?], where 


"i{R) — ||-Ri||p+ \\RkE\k-i\{.R 


Vllf- 




Under assumption 06.lip and with the stationarity of X, we derive that there 
exists C > 0 such that 


P ( sup 

\Si-l<p<Si 


Rn{p) -RRn{p) 


> ern] < N sup P(|i?„(r) - Ei?„(r)| > £r„) 

ri<r<r2 


< N- 

As N < c/ln and R = Pnb'/^, then for Vn = n 

n E sufficiently large. 


P ( sup 

Rn{p) - RRn{p) 

> £r A < 

\Si_l<p<Si 


) 

< 


KE 


-P 


rnb?,{y/n)P'^' 


KE 


-P 


^a(p+l)-d/2 ("jQg b(p+l) 


KE 


-P 


^a(p+l)-d/2-g2Qogn)^(P+^U'?l ' 

Eor a{p + 1) — d/2 — q 2 > 1 ei b{p + 1) + gi > 1, we get for any e: > 0 £ > 0 


n>l 


Si-l<p<Si 


Rn{p) - RRn{p) 


Considering these arguments the proofs of Theorem [3] are completed, it 
results from a direct application of the theorem of Borel-Cantelli and by 
Theorem [T] we have 

sup \RRn{r) — R{r)\ = 0{bf/) as n —)■ cxd. 

ri <r<r2 
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